ind_shnormMPEC <- function(expmeanval,expmu, oo, sharesum, marketForProducts) {

List<-ind(expmeanval, expmu, oo, sharesum)
denom1<-List$sum11[marketForProducts,]       # this expands the denominator
simShare<-(List$numer)*denom1                  # simulated shares for each draw
EstShare<-rowMeans(simShare) 		
EstShare<-as.matrix(EstShare)
simEstShare<-cbind(simShare,EstShare)
  
return(simEstShare)  

}
